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We argue that the analysis of the J-Q model, presented in Ref. and based on a field-theory 
description of coupled dimers, captures properly the strong quantum fluctuations tendencies, and 
the objections outlined in Ref. [I] are misplaced. 
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In a recent paper, Isaev, Ortiz, and Dukelsky 1 have 
questioned the interpretation of our results on the J-Q 
model 2 -. In addition, they have argued that their hierar- 
chical mean-field (HMF) method provides a more accu- 
rate description of the phase diagram of the model^. In 
our view the conclusions of Refs. OH are highly question- 
able, and we address the two relevant issues below. 

(1.) First, let us put our results (Ref. 0) in the correct 
perspective. In the notation of Ref. [J, let Jk be the 2- 
spin exchange and —K be the 4-spin one; both Jk and K 
are assumed positive. These are the quantities we have 
used in our work, Ref. [U (where we call Jk simply J). 
They are related to the couplings Q, J used in the original 
work of Sandviki via Q/J = (K/J K )/[1 - K/(2J K )]. 

We prefer the Jk,K units, because this way we can 
make the ratio of the 4- to the 2-spin term arbitrar- 
ily large, whereas, for technical reasons, in the origi- 
nal work of Sandvik this ratio has a maximum value of 
(Kj Jk) — 2, which corresponds to Q/J = oo. If (Kj Jk) 
is larger than 2, then this would mean J < 0, but still 
the 2-spin exchange in those units, (J+Q/2), is positive. 
Now, Sandvik's Monte Carlo (MC) result (confirmed also 
by later work^), gives the critical value (Q/J) c ~ 25, 
meaning (K/Jk) c ~ 1.85. Our most advanced calcula- 
tion finds the critical point at (K/Jk) c ~ 2.16, which 
we argue to be an improvement relative to our simple 
mean-field value, which is (K/ ' Jk)^ 1F % 1< 

Our main point is that while our mean-field result gives 
(K/J K )c IF £ 1 (i.e. Q/J small, less than 2), our im- 
provements, which take into account the strong quan- 
tum fluctuations, lead to {K/Jk)c around 2, meaning 
that \Q/J\ > 1. We call these "weak", and "strong" 
coupling regimes, respectively. In this sense our value 
(K/Jk) c = 2.16 is in "fairly good agreement with the 
MC" , i.e. the critical point is firmly in the strong- 
coupling regime. 

The authors of Ref. [l| object to the fact that our 
(K/Jk) c — 2.16 corresponds to {Q/J) c = —27, i.e. 
J < 0, which appears very far from the Monte Carlo: 
(Q/J) c = +25, i.e. (K/J K )c = 1-85. However this a 
very misleading and incorrect way to look at the results 



(that's why we prefer the K, Jk units). 

Since the transition is at large \Q/J\, it is irrelevant 
that J < 0. Indeed, the ratio of the 4- to the 2-spin term 
at our critical point, (K/Jk) c — 2.16, corresponds to the 
4-spin term being Q, and the 2-spin term being {—\J\ + 
Q/2), with Q/\J\ = 27. Sandvik's result, (K/J K ) C = 
1.85, means that the 4-spin term is Q, and the 2-spin 
term is (J + Q/2), with Q/J = 25. 

In both cases \Q/J\ 3> 1, so the sign of J is irrelevant. 
This is what we call the strong-coupling regime, when 
Kj Jk is around 2, and we describe the agreement be- 
tween {K/J K )c = 1-85 and {K/J K ) C = 2.16 as "fairly 
good". The two critical points are fairly close, when 
things are put in the right context. We have certainly not 
achieved a perfect agreement with the Monte Carlo, but 
we have found the correct quantum fluctuations trend. 

Of course, as pointed out in Ref. lj, the point 
[Kj Jk)c = 2.16 is outside the range explored in the work 
of Sandvik^, but the physics is expected to be the same, 
i.e. one simply penetrates deeper into the quantum dis- 
ordered, gapped phase. 

(2.) We also emphasize that our results were ob- 
tained under the assumption that the ground state is 
of the columnar dimer type, as argued in Ref. i. In 
our work we found that this ground state is stable for 
K/Jk > {K/Jk)c, but we did not compare with other 
possible ground states, as it is quite difficult to compare 
ground state energies reliably^, especially when the sys- 
tem exhibits strong fluctuations. 

A severe problem of the HMF approach^ is that the 
critical point is at (Q/J) c ~ 2, i.e. (K/Jk) c ~ 1, which 
is far off the Monte Carlo result. In fact the HMF crit- 
ical point location is close to the one we found in our 
"simple" dimer mean-field framework 2 , which we ruled 
out as unreliable. Thus the hierarchical plaquette mean- 
field seems to place the critical point firmly in the weak- 
coupling (small Q/J) region, in disagreement with the 
Monte Carlo, and consequently it is highly unlikely that 
the HMF takes properly into account the strong quantum 
fluctuations present in the J-Q model. 
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